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The Principle of Maximum Conformality (PMC) eliminates QCD renormalization scale-setting 
uncertainties using fundamental renormalization group methods. The resulting scale-fixed pQCD 
predictions are independent of the choice of renormalization scheme and show rapid convergence. 

The coefficients of the scale-Hxed couplings are identical to the corresponding conformal series with 
zero /3-function. Two all-orders methods for systematically implementing the PMC-scale setting 
procedure for existing high order calculations are discussed in this article. One implementation is 
based on the PMC-BLM correspondence (PMC-I); the other, more recent, method (PMC-II) uses the 
72.i-scheme, a systematic generalization of the minimal subtraction renormalization scheme. Both 
approaches satisfy all of the principles of the renormalization group and lead to scale-fixed and 
scheme-independent predictions at each finite order. In this work, we show that PMC-f and PMC- 
II scale-setting methods are in practice equivalent to each other. We illustrate this equivalence for 
the four-loop calculations of the annihilation ratio R„+f,- and the Higgs partial width T{II —>■ bb). 

Both methods lead to the same resummed (‘conformal’) series up to all orders. The small scale 
differences between the two approaches are reduced as additional renormalization group {/3i}-terms 
in the pQCD expansion are taken into account. We also show that special degeneracy relations, 
which underly the equivalence of the two PMC approaches and the resulting conformal features of 
the pQCD series, are in fact general properties of non-Abelian gauge theory. 

PACS numbers: 12.38.Bx, 12.38.Aw, 11.15.Bt 


A primary problem for perturbative QCD (pQCD) 
analyses is how to systematically set the renormaliza¬ 
tion scale of the QCD running coupling to achieve pre¬ 
cise fixed-order predictions for physical observables [1]. 
A valid prediction at any finite order should be indepen¬ 
dent of the choice of the renormalization scheme, since 
the choice of the scheme is a theoretical convention. The 
“Principle of Maximum Conformality” (PMC) [2-5] elim¬ 
inates the renormalization scheme- and scale- ambiguities 
order by order in perturbation theory, consistent with 
fundamental renormalization group methods [6-9]; it is 
also the principle underlying the well-known Brodsky- 
Lepage-Mackenzie (BLM) approach [10]. In the Abelian 
limit, the PMC is the standard scale-setting method used 
for precision tests of quantum electrodynamics (QED). 
The elimination of the renormalization scale ambigu¬ 
ity removes an unnecessary systematic error for pQCD 
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predictions, thus providing scheme-independent precision 
tests of the Standard Model and improving the sensitivity 
to new physics. 

The renormalization group equation (RG-equation) 
determines the running of the gauge coupling from the 
analytic properties of the /3-function: 


l3{as{R)) = ( 1 ) 

^ i=0 

where a perturbative expansion of the /3-function in terms 
of the coupling Us = asj47r is assumed. The expres¬ 
sions for Pq, ..., Ps in the modified minimal-subtraction 
scheme (MS-scheme) can be found in Refs. [11-15]. 
One can then use a Taylor-expansion to derive a scale- 
displacement relation for the running coupling at two dif¬ 
ferent scales /xi and ^ 2 ' 
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The RG-equation and the asymptotically expanded /3- 
function in (1) can be used to recursively establish a per- 
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turbative /3-patterii at each order: 
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order [2-5]. Both implementations have been success¬ 
fully used for making scale-fixed and rapidly converging 
predictions for a number of high-order high-energy pro¬ 
cesses [17-26]. We shall show that the two PMC meth¬ 
ods are, in practice, equivalent. Any small residual scale 
difference between the two approaches is systematically 
reduced as more RG {/3i}-terms in the perturbative QCD 
expansion are known. Both implementations of the PMC 
allow one to determine the PMC scales of a process order- 
by-order in pQCD. 

In general, a pQCD prediction for a physical observable 
can be written as an expansion ^ 


The PMC utilizes this perturbative /3-pattern to system¬ 
atically set the scales of the running coupling at each 
order in a pQCD expansion; the coefficients of the result¬ 
ing series thus match the coefficients of the corresponding 
conformal theory with /3 = 0. Thus the divergent a"/3”n! 
“renormalon” contributions are absorbed into the scale- 
fixed running couplings. The pQCD convergence of the 
resummed series is generally improved due to the elimi¬ 
nation of those renormalon terms. This is the same prin¬ 
ciple used in QED where all /3-terms resulting from the 
vacuum polarization corrections to the photon propaga¬ 
tor are absorbed into the scale of the running coupling. 
As in QED, the scales are physical in the sense that they 
reflect the virtuality of the gluon propagators at a given 
order, as well as setting the effective number nj of active 
flavors. The resulting resummed pQCD expression thus 
determines the relevant “physical” scales for any physi¬ 
cal observable, thereby providing a solution to the renor¬ 
malization scale-setting problem. There can be a small 
residual scale-uncertainty in the final expression due to 
the truncation of the /3-function; however, these residual 
uncertainties are highly suppressed even for lower-order 
predictions. 

The scale-displacement relation (2) provides the sim¬ 
plest example of the PMC: If one resums all-orders {Pi}- 
terms in the right-hand-side of Eq.(3), one (trivially) ob¬ 
tains as(/J. 2 ) = as(/r 2 )lpMC at any fixed order, indepen¬ 
dent of /ii and the choice of scheme; this agrees with the 
reflexivity property of the renormalization group [16]. 

In this paper we will show how to systematically imple¬ 
ment the PMC scale-setting if one starts with an existing 
high-order pQCD calculations for a physical observable. 
Such calculations which are available in the literature 
are usually performed in a conventional renormalization 
scheme such as the MS-scheme, assuming a single initial 
scale /r. Two all-orders approaches for implementing the 
PMC procedure for existing high-order calculations have 
been suggested. One is based on the PMC-BLM corre¬ 
spondence (PMC-I) [2]; the other method (PMC-II) [4, 5] 
utilizes the pattern of P terms illuminated by the TZg- 
scheme, which was introduced in Ref. [4] as a general¬ 
ization of the MS-like renormalization schemes, thus en¬ 
abling one to obtain nontrivial information on the pQCD 
series. Both approaches satisfy all of the principles of 
the renormalization group [16], and they lead to scale- 
fixed and scheme-independent predictions at any hnite 
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where the symbol “...” means even higher-order contri¬ 
butions, /X is the initial scale. The n/-terms count the 
number of active flavors which arise from light-quark loop 
contributions to the /3-function. All of the explicitly writ¬ 
ten n/-terms thus pertain to the RG {/3i}-terms, and - 
as in the BLM procedure - they provide a guide to set¬ 
ting the renormalization scales. It is important to note 
that the coefficients may also contain n/-terms which 
are ultraviolet finite and are not associated with the /3- 
function. Eor example, terms arising from the ultravio¬ 
let finite quark-loop contributions to the three-gluon and 
four-gluon vertices will be kept unchanged during PMC 
scale-setting. 

After applying PMC-I and PMC-II, all of the n/-terms 
which are governed by the RG-equation will be resummed 
into the running coupling; the pQCD series in Eq. (4) 
then changes to the resummed ‘conformal’ series: 

m 

P = .) + ..., (5) 

i=l 

where M=I or II for PMC-I or PMC-II, respectively. The 
coefficients are the /3-independent ‘conformal’ coeffi¬ 
cients and Qm.i are the PMC scales, where i = 1 stands 
for the leading-order (LO) one, i = 2 stands for the next- 
to-leading (NLO) on, etc. The resulting coefficients are 
conformal in the sense that they are UV finite, scheme- 
independent, and do not depend on the resummation 
schemes, PMC-I, PMC-II, or any other implementation 
following PMC. This will be shown explicitly for PMC-I 
and PMC-II. 

Let us describe the different implementations of PMC 
in more detail: 

PMC-I allows one to obtain the correct PMC scales us¬ 
ing a step-by-step method without first transforming the 


^ We do not consider quark mass renormalization in this paper and 
will assume n > 1 in our discussions. For the case of n = 0, i.e. 
when the tree-level result does not involve strong interactions, 
one can implement the PMC from the second a^-term and all 
the following formulas apply. 
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n/-terms into the {/3i}-terms. This procedure is based 
on the observation that one can rearrange all the Feyn¬ 
man diagrams of a process in form of a cascade; i.e., the 
“new” terms emerging at each order can be equivalently 
regarded as a one-loop correction to all the “old” lower- 
order terms. All of the n/-terms can then be absorbed 
into the running coupling following the basic /3-pattern in 
the scale-displacement formula, i-e. Eq.(3). The resulting 
PMC-I scales are themselves expressed as perturbative 
expansions with the same /3-pattern of Eq.(3). More ex¬ 
plicitly, the PMC-I scales can be derived in the following 
way: The LO PMC-I scale Qpi is obtained by eliminat¬ 
ing all the n/-terms with the highest power at each order, 
and at this step, the coefficients of the lower-power n/- 
terms are changed simultaneously to ensure that the cor¬ 
rect LO cTs-running is obtained; the NLO PMC scale Qp 2 
is obtained by eliminating the ny-terms of one less power 
in the new series obtaining a third series with less n/- 
terms; and so on until all n/-terms are eliminated. The 
step-by-step coefficients for the n/-series can be found 
in Ref. [2]. After performing these order-by-order scale 
shifts /i —>■ Qip, Qi,i —t (5p2i Qi ,2 ( 5 i, 3 , • • •, one elim¬ 
inates all the n/-terms associated with the ag-running 
and derive the conformal series. 

It is noted that the PMC scale-setting can be automat¬ 
ically implemented in a higher order pQCD calculation if 
one uses the 72.5-scheme. The usual subtraction constant 
IndTT — 7 e is generalized to In dvr — 7 ^ — i5. The resulting 
dependence on the extra constant 5 flags all of the terms 
in the pQCD series proportional to the /3-function [4, 5]. 
Thus, in contrast to PMC-I, PMC-II first transforms the 
n/-series at each order into the specific ,5-pattern dic¬ 


tated by the T^^-scheme; the resulting /3-pattern leads 
directly to the PMC scales and the conformal series. In 
this sense, PMC-II is a theoretical improvement of PMC- 
1, since PMC-I only determines the overall effective scales 
without determining the perturbative terms which lead 
to those values. Due to the fact that the running of ttg 
at each order has its own {,5i }-series as governed by the 
RG-equation, the /3-pattern for the pQCD series at each 
order is a superposition of all of the {/3i}-terms which 
govern the evolution of the lower-order contributions 
at this particular order. The resulting /3-pattern is then 
in general different from the /3-pattern of Eq.(3). 

PMC-II suggestes that the coefficients of the {/3i}- 
terms in the /3-pattern can be fixed by requiring a “de¬ 
generacy relation” among different {,5i}-terms at differ¬ 
ent orders; the result resembles a skeleton-like expan¬ 
sion [4, 5]. By resumming the {/3i}-series according to 
this expansion, one also correctly reproduces the Abelian 
A’c —>■ 0 limit of the observables [27]. Thus one can simul¬ 
taneously determine the PMC scales Qny at all orders 
from their initial values /i; i.e. p, —?► Qii,i, t Qii, 2 , 
Ai —t Qii.s, •••, by resumming the {/3i}-terms into the 
running couplings in the skeleton-like form. 

The degeneracy relations introduced by PMC-II were 
originally obtained by studying the pQCD series in the 
725 -scheme [4]. Let us show that theses relations are re¬ 
quired by the conformality of the final series. Using the 
RG-equation and the scale displacement relation, we can 
write down the most general /3-pattern for the pQCD ap- 
proximant of a physical observable. For example, Eq.(4) 
can be rewritten as 


P = + (^"0 + -f + {n + l)/3orf^i + ^ 
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By applying the PMC-II procedures [4, 5], we obtain 

P = rl^<iQiip) + rl, • a:+\Qna) + + n - rl,] 

+ {nl32 Kp* - r^h] -f (n -f I)/3i -h [f"; - ^ a^s""^ip) + ■■■ ■ 


( 6 ) 


(7) 


This step does not require the degeneracy relations. How¬ 
ever, in order to ensure that p is conformal; i.e., to not 
contain any {/3j}-terms, we get the required degeneracy 
relations: rf l* = ^ etc.. 

Those degeneracy relations also ensure the elimination 
of all the remaining as{p)-teTm.s in the pQCD series, we 


then obtain the required scale-fixed PMC prediction. 

Alternatively, if one implements PMC-I and requires 
the uniqueness of the conformal coefficients for the two 
methods, we again find the degeneracy relations. To show 
this, we transform the general {/3i}-series ( 6 ) back to n/- 
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series, and by further applying PMC-I, one finds 

„i _ ^ii 
' 1,0 “ ' 1 , 0 ! 

„i _ 

' 2,0 “ ' 2 , 0 ! 
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^i,o = riio-icl(n(151riii- 228riit + 77r"r)CF 

-21(n + l)(r^^i - + ^n(41r^^i + 120r^^i 


- Y^ACF{n{7rl, - 


Thus to ensure the conformality of the PMC-I and the 
PMC-II final expressions; i.e. i"Jq = rpg, we are immedi¬ 
ately led to the degeneracy relations. Following the same 
procedures, we can demonstrate the equivalence of con¬ 
formal series up to any order. This equivalence can be 
explained by the fact that the scale-displacement relation 
( 2 ) acts only on the purely non-conformal {/3i}-series, 
and PMC-I and PMC-II only differ in eliminating the 
n/-terms - either by using the RG /3-pattern directly, or 
by using the super-positioned RG /3-pattern. Thus, if one 
transforms the PMC-I prediction to the one of PMC-II, 
or vice versa, the conformal coefficients are not altered. 
More explicitly, by using the /3-function to the highest 
known order; i.e., four-loops, the conformal coefficients 
for any semi-simple Lie gauge group with n/-fermions 
and Nc-co\ots are: 


I.II 
'’1,0 

I.II 
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16T|c3,o + 44CaTfC3,i + I2ICIC3.2] , (10) 

^[2CaT^{-287CI + 1208CaCf + 

924C|)c2,i - ^8CaTI{7Ca + 11Cf)c3,i - 
264C'irF {7Ca + HCf) C3,2 + 64T|c4,o + 
l7QCATlciA + 484C'iTFC4,2 + 1331C'ic4,3], 

nil 


of the PMC-I and PMC-II scales by = InQpi/Qipi, 
and the first three ones for a four-loop prediction read: 


Ai = 
A 2 = 

A 3 = 
where 


3/3i/ 
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/ = 6n^(n-I-l)ci qC 4_3-I-(2n^-I-8n^-I-13n-I- 7 ) C2_i 
—6n (n^ -I- 3n -I- 3 ) ci,oC3,2C2,i, 
g = 4n(n-f l)(n-I-2)ci^oC2,i(C'^(84c2,irF — 12103 ^ 2 ) 
+44CaTf(3c2,iCf - 03 , 1 ) - 16c3,oT|). 


As indicated by Eqs.(12,13,14), the LO logarithmic scale 
difference Ai starts at the order a1, which changes to 
order for the NLO A 2 , and to order a° for the NNLO 
A 3 . The value of A^ can be qualitatively understood 
through the scale-displacement equation (3), i.e., 

a\{Qii,i) - <(Qi,i) = 2i/3oAia*+i((5n,i) + 0{a\+‘^). 

The leading term of A^ may be of order a°. Since the 
NLO n/-term is uniquely fixed, the order a° of InQi 1 
and InQii^i must be equal. This explains why PMC- 
1, PMC-II, and also BLM, are exactly the same at the 
NLO level [29]. We further note that the order al of 
InQi.i and InQn,!, and the order oflnQi 2 and lnQii .2 
are also equal. This shows a non-trivial equivalence of 
the two scale-setting methods at the non-conformal level, 
which means that PMC-I and PMC-II are the same at the 
NNLO level. Moreover, if we know additional higher-loop 
contributions, we can achieve more precise and closer 
PMC-I and PMC-II scales, and thus obtain smaller A^. 



Al 

A 2 

A 3 

Re+e- (Q = 

31.6 GeV) 

-0.0043 

+0.0973 

+1.9389 

pNS/^ ^ 

hadrons) 

-0.0030 

-0.0826 

+2.0432 

r(T(is) 

-A e+e") 

+0.0353 

+0.1047 

+0.0816 

r{H 

bb) 

+ 0.0001 

-0.0014 

-0.0018 


where Ca, Cf and Tp are quadratic Casimir invari¬ 
ants [28]. For a SU{Nc)-co\oy group, we have Ca = Nc, 
Cf = {N^ - l)/2Ac and Tp = 1/2. 

Since the non-conformal {/3i}-terms are eliminated in 
different ways, the PMC-I and PMC-II scales can in prin¬ 
ciple be different. Since they are both based on the RG- 
equation, the accuracy of the PMC scales depend heavily 
on how well we know the {/3i}-terms of the process. One 
observes that the PMC-I and PMC-II scales are them¬ 
selves expressed as perturbative series, and their logarith¬ 
mic differences will be suppressed by specific powers of 
Us- To quantify this, we define the logarithmic difference 


TABLE I: The logarithmic scale difference Ai for i?e+e- {Q = 
31.6GeV), r^^{Z -A hadrons), F (T(IS')e+e"), and 
T{H -A bb) np to four-loop QCD corrections, where NS stands 
for the non-singlet contribution. 

Table I shows several four-loop examples of how Ai 
changes with the increment of the loop corrections. The 
four-loop expressions using conventional scale setting are 
adopted from Refs. [30-37]. The PMC predictions for 
those channels can be found in Refs. [23-26]. In the case 
of 77 —>• bb, the computed scale differences are very small 
at any order. For example, the largest difference for the 
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case oi H ^ bb is found between and Qii^a, which 
is less than 0.2%. The logarithmic scale differences for 
Rie + e—) and —>■ hadrons) are somewhat larger; 

their magnitudes follow the trend |Ai| ^ IA 2 I IA 3 I, 
indicating that the PMC-I and PMC-II scale differences 
quickly diminish as we include more {/?i}-terms to fix the 
PMC scales. In the case of P (T(IS') —>• e+e“), the log¬ 
arithmic scale differences |Ai| ^ IA 2 I IA 3 I show that 
the scale differences between Qi 3 and Qii.s are acciden¬ 
tally small even with less /3-term information than the 
case of t^i/ii. 2 - 

In the following, we will consider two observables 
Re+e-{Q = 3I.6GeV) and T{H —>• bb) in detail to illus¬ 
trate the differences and common features of the PMC-I 
and PMC-II predictions. 

We rewrite the e“'"e“ annihilation i?-ratio and the 
H ^ bb decay width as 

Re+e-iQ) = 3^e2(l + i?))(Q,^)), (15) 

Q 

r{H^bb) = (l + J?,f(g,/r)), (16) 

where = I]r=o ^he 

scale fj, is arbitrary. Here Q stands for a typical mo¬ 
mentum flow of the process. In order to compare with 
data, we will take Q = 31.6 GeV for the i?-ratio and 
Q = Mh for H —>■ bb. The coefficients with 

its explicit ny-dependence up to four-loop level can be 
found in Refs. [30, 31]. We adopt four-loop as-running 
and fix the QCD parameter by using as{Mz) = 

0.1185 ±0.0006 [38]. 



LO 

NLO 

N^LO 

N^LO 

total 

PMC-I 

0.04294 

0.00340 

-0.00002 

-0.00001 

0.04631 

PMC-II 

0.04290 

0.00352 

-0.00004 

-0.00002 

0.04636 

Conv. 

0.04499 

0.00285 

-0.00117 

-0.00033 

0.04634 


TABLE II: The contributions of each of the contributions 
(LO, NLO, N^LO and N^LO) to the four-loop pQCD ap¬ 
proximate i? 3 . 



LO 

NLO N^LO N^LO 

total 

PMC-I 

0.2268 

0.0249 -0.0091 -0.0012 

0.2414 

PMC-II 

0.2268 

0.0249 -0.0094 -0.0012 

0.2411 

Conv. 

0.2037 

0.0377 -b0.0019 -0.0014 

0.2419 


TABLE III: The contributions of each of the contributions 
(LO, NLO, N^LO and N^LO) to the four-loop pQCD ap¬ 
proximate . 

Tables II and III show the contributions of each con- 

e/H 

tribution up to the four-loop pQGD term R^' . The 

pQCD predictions using conventional scale setting with 
fi = Q are also presented for comparison. We have set 
the N^LO scale Qp 4 (Qii, 4 ) to the be highest order de¬ 
termined scale Qi ,3 (Qips), following the prescription of 
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FIG. 1: Results for f?® and up to (n-l-l)-loop QCD correc¬ 
tions together with their approximate errors [max 

under the PMC-I and PMC-II methods, where is the 

(n + l)th-order conformal coefficient and n = (1, 2, 3). 


PMC-I. Tables II and III show that the pQCD conver¬ 
gence for the PMC-I and PMC-II are very similar, which 
are, as required, better than the conventional one due 
to the elimination of divergent renormalon terms. The 
pQCD series for H ^ bb is almost the same. There 
are slight differences for R{e~^e~) due to the logarith¬ 
mic scale differences listed in Table I; however, these are 
greatly suppressed by the magnitudes of the conformal 
coefficients and the ag-powers. 

FIG. 1 presents the results for Rn^ up to (n ± I)- 
loop QGD corrections, together with their estimated er¬ 
rors ±lCn^^a"+^lMAX using both the PMC-I and PMC- 
II approaches. Here, Cn^ are (n ± l)th-order conformal 
coefficients, and the subscript “MAX” means the maxi¬ 
mum value of lCn^^a"+^j obtained by varying the scale 
/i within the usual region of [Q/2,2Q]. This error esti¬ 
mate is natural for the PMC, since after applying the 
PMC-I or PMC-II, the pQCD convergence is ensured 
and the only uncertainty is from the last term of the 
pQCD series due to the unfixed PMC scale at this par- 
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ticular order. When additional pQCD loop corrections 
are taken into consideration, one obtains a weaker scale 
dependence. This agrees with the conventional wisdom 
that as one incorporates a higher order calculation, one 
can obtain an increasingly reliable scale-invariant esti¬ 
mate. The PMC-I and PMC-II predictions are very close 
in magnitude; their values quickly approach convergence, 
indicating that a low-order calculation could be sufficient 
to achieve an accurate pQCD prediction. 

Let us end with a final comment on the general¬ 
ity of the degeneracy relations found using the TZs- 
scheme [4] . The degeneracy relations actually hold under 
any scheme, which was pointed out in Ref. [5]. To explain 
this, we adopt the effective charge method introduced by 
Grunberg [39]. Any effective charge a a of a physical ob¬ 
servable A can be expanded over the running coupling 
an of the 7?.5-scheme as 

aA{^J) = + [r^o + Poril]a\{^l) + [r^o + 

/3ir^i -h 2/3or^i -h /3o^^2]a|i(Ai) + 

Ko + + 2/3ir^i -h + 

+ 3/3oT’^2 + 0o'rt‘i]an{lJ-f H-(A7) 

where the universality of /3o and /3i is used, and the 
explicit scheme-dependence of (32 is expressed. The (3- 
function of the A-scheme is related to the one of the Rg- 
scheme through the identity, (3^{a a) = If 

one analyzes another effective charge ob in the same way, 
one finds that a a can also be expanded over as, through 
scheme-transformations mediated by the R^-scheme 

aAip) = + {r2,o + + 

+ 2/3of;yf + + 

(rto + + 2/3ir^f + ^l3ol3ir^^ + 

+ ■ ■ -( 18 ) 


This shows that the degeneracy relations still hold, even if 
scheme B is not an MS-like scheme. The coefficients 
up to four-loop level can be found in Ref. [5] Among 
them the conformal ones are 
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= 
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2,0 


2 r. 


(19) 

( 20 ) 

3'!o)^2^o(;21) 


which are purely expressed in terms of the conformal co¬ 
efficients in the T^-^-scheme. This self-consistency 

shows that the applicability of the PMC is ensured and 
its predictions are scheme-independent. 

In conclusion, we have shown that the two all-orders 
PMC approaches are equivalent to each other at the level 
of conformality and are thus equally viable PMC proce¬ 
dures. PMC-I implementation is a direct extension of 
the BLM approach, whereas PMC-II provides additional 
theoretical improvements; in addition, it can be read¬ 
ily automatized using the i? 5 -scheme. By construction, 
both the PMC-I and PMC-II satisfy all of the principles 
of the renormalization group, thus providing scale-fixed 
and scheme-independent predictions at any fixed order. 
Those two implementations of PMC differ, however, at 
the non-conformal level, by predicting slightly different 
RG scales of the running coupling. This difference arises 
due to different ways of resumming the non-conformal 
terms, but this difference decreases rapidly when addi¬ 
tional loop corrections are included. 


The key step of PMC-II is to use the pattern gener¬ 
ated by the RG-equation and its degeneracy relations to 
identify which terms in the pQCD series are associated 
with the QCD /3-function and which terms remain 
in the /3 = 0 conformal limit. The /3-terms are then 
systematically absorbed by shifting the scale of the 
running coupling at each order, thus providing the PMC 
scheme-independent prediction. The recursive patterns 
and degeneracy relations between the /3-terms at each 
order are essential for carrying out this procedure. 
The implementation of PMC-II illuminates how the 
renormalization scheme and initial scale dependence 
are eliminated at each order. These advantages shows 
the PMC-II is theoretically robust and is the preferred 
method for practical implementations of PMC. 


Acknowledgements: This work was supported in part 
by the Natural Science Foundation of China under 
Grant No.11275280, the Fundamental Research Funds for 
the Central Universities under Grant No.CDJZR305513, 
and the Department of Energy Contract No.DE-AC02- 
76SF00515. SLAC-PUB-16287. Nordita-2015-54. 


[1] X. G. Wu, S. J. Brodsky and M. Mojaza, Prog. Part. 
Nucl. Phys. 72 , 44 (2013). 

[2] S. J. Brodsky and X. G. Wu, Phys. Rev. D 85, 034038 
(2012) [Erratum-ibid. D 86, 079903 (2012)]. 

[3] S. J. Brodsky and L. Di Giustino, Phys. Rev. D 86, 
085026 (2012). 


[4] M. Mojaza, S. J. Brodsky and X. G. Wu, Phys. Rev. 
Lett. 110 , 192001 (2013). 

[5] S. J. Brodsky, M. Mojaza and X. G. Wu, Phys. Rev. D 
89, no. 1, 014027 (2014). 

[6] E. G. G. Stueckelberg and A. Petermann, Helv. Phys. 
Acta 26 , 499 (1953). 



7 


[7] A. Peterman, Phys. Kept. 53, 157 (1979). 

[8] M. Gell-Mann and F. E. Low, Phys. Rev. 95, 1300 (1954). 

[9] N.N. Bogoliubov and D.V. Shirkov, Dok. Akad. Nauk 
SSSR 103, 391 (1955). 

[101 S. J. Brodsky, G. P. Lepage and P. B. Mackenzie, Phys. 
Rev. D 28, 228 (1983). 

[11] H. D. Politzer, Phys. Rev. Lett. 30, 1346 (1973). 

[12] O. V. Tarasov, A. A. Vladimirov and A. Y. Zharkov, 
Phys. Lett. B 93, 429 (1980). 

[13] T. van Ritbergen, J. A. M. Vermaseren and S. A. Larin, 
Phys. Lett. B 400, 379 (1997). 

[14] K. G. Ghetyrkin, Nucl. Phys. B 710, 499 (2005). 

[15] M. Czakon, Nucl. Phys. B 710, 485 (2005). 

[16] S. J. Brodsky and X. G. Wu, Phys. Rev. D 86, 054018 

( 2012 ). 

[17] S. J. Brodsky and X. G. Wu, Phys. Rev. D 86, 014021 
(2012) [Erratum-ibid. D 87, 099902 (2013)]. 

[18] S. J. Brodsky and X. G. Wu, Phys. Rev. D 85, 114040 

( 2012 ). 

[19] S. J. Brodsky and X. G. Wu, Phys. Rev. Lett. 109, 
042002 (2012). 

[20] S. Q. Wang, X. G. Wu, Z. G. Si and S. J. Brodsky, Phys. 
Rev. D 90, 114034 (2014). 

[21] S. Q. Wang, X. G. Wu, J. M. Shen, H. Y. Han and Y. Ma, 
Phys. Rev. D 89, 116001 (2014). 

[22] S. Q. Wang, X. G. Wu, X. G. Zheng, G. Ghen and 
J. M. Shen, J. Phys. G 41, 075010 (2014). 

[23] S. Q. Wang, X. G. Wu and S. J. Brodsky, Phys. Rev. D 
90, 037503 (2014). 

[24] S. Q. Wang, X. G. Wu, X. C. Zheng, J. M. Shen and 


Q. L. Zhang, Eur. Phys. J. G 74, no. 4, 2825 (2014). 

[25] X. G. Wu, Y. Ma, S. Q. Wang, H. B. Fu, H. H. Ma, 
S. J. Brodsky and M. Mojaza, arXiv: 1405.3196 [hep-ph]. 

[26] J. M. Shen, X. G. Wu, H. H. Ma, H. Y. Bi and 
S. Q. Wang, arXiv: 1501.04688 [hep-ph]. 

[27] S. J. Brodsky and P. Huet, Phys. Lett. B 417, 145 (1998). 

[28] T. van Ritbergen, A. N. Schellekens and J. A. M. Ver¬ 
maseren, Int. J. Mod. Phys. A 14, 41 (1999). 

[29] H. H. Ma, X. G. Wu, Y. Ma, S. J. Brodsky and M. Mo¬ 
jaza, arXiv: 1504.01260 [hep-ph]. 

[30] P. A. Baikov, K. G. Ghetyrkin and J. H. Kiihn, Phys. 
Rev. Lett. 96, 012003 (2006). 

[31] P. A. Baikov, K. G. Ghetyrkin and J. H. Kiihn, Phys. 
Rev. Lett. 101, 012002 (2008). 

[32] M. Beneke, Y. Kiyo and K. Schuller, Phys. Lett. B 658, 
222 (2008) 

[33] M. Beneke, Y. Kiyo and K. Schuller, Nucl. Phys. B 714, 
67 (2005) 

[34] A. A. Benin, V. A. Smirnov and M. Steinhauser, Nucl. 
Phys. B 716, 303 (2005) 

[35] M. Beneke, Y. Kiyo and A. A. Benin, Phys. Lett. B 653, 
53 (2007) 

[36] M. Beneke and Y. Kiyo, Phys. Lett. B 668, 143 (2008) 

[37] P. A. Baikov, K. G. Ghetyrkin, J. H. Kiihn and J. Rit- 
tinger, Phys. Rev. Lett. 108, 222003 (2012) 

[38] K. A. Olive et al. [Particle Data Group Gollaboration], 
Chin. Phys. C 38, 090001 (2014). 

[39] G. Grunberg, Phys. Rev. D 29, 2315 (1984). 


